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ABSTRACT 

A  statistical  analysis  of  the  "Second  Catalog  of 
Emission- Line  Stars  of  the  Orion  Population"  was  undertaken. 
After  the  reduction  of  the  Catalog  to  numeric  data,  statis- 
tics and  frequency  distributions  were  systematically  searched 
for  any  physically  significant  characteristics.  Those  two-way 
frequency  distributions  with  possible  physical  significance 
were  analyzed. 

Distance  moduli  for  228  stars  were  tabulated. 
Hertzsprung-Russell  diagrams  were  plotted.   The  T  Tauri  stars, 
so  plotted,  lie  between  the  main  sequence  and  giant  branch. 
A  lack  of  correlation  between  the  large  scale  variability 
and  any  other  set  of  observed  characteristics  was  discovered. 
Present  theories  of  stellar  evolution  do  not  explain  this 
lack  of  correlation  nor  the  cause  of  the  variability.   It 
is  argued,  that  on  the  basis  of  current  data,  no  conclusion 
can  be  reached  as  to  whether  T  Tauri  stars  are  on  convective 
or  radiative  equilibrium  tracks  approaching  the  main  sequence. 
A  luminosity  function  for  the  Orion  Population  is  presented 
and  is  found  to  be  significantly  different  than  the  solar 
neighborhood  luminosity  function. 
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I.   INTRODUCTION 

The  dynamics  and  evolution  of  stars  in  hydrostatic 
equilibrium  fusion  cores  is  well  understood.   There  is 
relatively  good  agreement  between  theory  and  observations 
on  stellar  age,  composition,  luminosity  and  effective  tem- 
perature, and  stability  of  these  objects.   Many  observations 
support  these  theories  that  predict:  shorter  lifetimes  for 
more  massive  stars;  evolution  of  giant  and  supergiant  stars; 
production  of  planetary  nebulae,  novae,  and  super  novae;  and 
collapse  to  white  dwarfs  and  neutron  stars.   It  is  possible 
that  even  the  final  stage  of  stellar  evolution,  (for  some 
stars),  the  black  hole,  has  been  indirectly  observed. 

No  such  set  of  matching  observations  and  theories  exist 
for  newly  formed  stars.   Several  theories  have  been  suggested 
which  predict  the  details  of  collapse  for  dense  clouds  of 
interstellar  gas  and  dust,  under  highly  idealized  conditions. 
But  observational  support  for  many  of  the  gross  physical 
details  of  these  schemes  has  not  been  found.   One  prediction 
of  the  theories  of  early  stellar  evolution  is  the  appearance 
of  luminous  young  stars.   The  stars  at  the  time  of  the  appear 
ance  are  predicted  to  have  luminosities  greater  than  main 
sequence  stars  of  the  same  effective  temperature.   If  stars, 
whose  pre-main  sequence  nature  can  be  determined  are  studied, 
then  their  observed  properties  can  be  used  to  distinguish 
between,  or  even  challenge  the  basic  features  of  existing 
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theories  of  early  stellar  evolution.   The  T  Tauri  stars 
seem  to  be  such  a  group,  by  reason  of  their  Hcrtzsprung- 
Russell  diagram  position,  spectral  peculiarities,  and  other 
unusual  properties. 

The  first  recognition  of  the  homogeneous  group  of  irreg- 
ular variable  (in  apparent  magnitude)  stars  was  made  by  Joy 
[1945].   Eleven  stars  were  included  in  that  group,  including 
the  prototype  star  T  Tauri,  for  which  the  group  is  named. 
Although  no  two  of  the  T  Tauri  stars  are  exactly  alike,  Joy 
noted  four  characteristics  that  distinguish  T  Tauri  variables 
from  the  multitude  of  other  variable  stars.   These  are:  "(1) 
rapid  irregular  light- variations  of  about  three  magnitudes; 
(2)  spectral  type  F5-G5,  with  emission  lines  resembling  those 
of  the  solar  chromosphere,  particularly  in  the  great  strength 
of  H  and  K  of  calcium;  (3)  low  luminosity;  (4)  association 
with  dark  or  bright  nebulosity."   In  intervening  years, 
these  stars  have  come  to  be  accepted  as  pre-main  sequence 
obj  ects . 

Objective  prism  surveys  since  then  have  increased  the 
number  of  known  T  Tauri  stars  to  more  than  three  hundred, 
including  many  stars  which  because  of  their  association  with 
nebulosity,  emission- 1 ine  spectra,  and  irregular  light  vari- 
ations are  T  Tauri-like  even  though  they  are  not  strictly 
T  Tauri  stars  by  Joy's  criteria.   Herbig  [1962]  suggested 
that  this  overall  set  of  emission  and  irregular  variable 
stars  be  referred  to  as  the  "Orion  Population."   By  1972, 
Herbig  and  Rao  [1972]  were  able  to  catalog  323  stars  in  the 
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"Orion  Population,"  most  of  which  are  T  Tauri  stars.   Al- 
though the  information  is  incomplete  for  some  stars,  the 
Catalog,  as  the  Herbig-Rao,  1972,  paper  will  hereafter  be 
called,  includes  data  for  the  categories  listed  in  Table  I. 
Concurrently  in  the  table,  the  Fortran  variable  name  assigned 
to  the  particular  characteristic  is  listed.   These  Fortran 
names  are  sometimes  used  in  later  discussions  of  the  Catalog 
data.   The  program  undertaken  was  to  calculate  the  distribu- 
tions and  relationships  among  the  categories  in  the  Catalog. 
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II.   BACKGROUND  OF  THE  T  TAURI  PHENOMENA 

A  deeper  discussion  of  the  properties  of  member  stars 
of  the  Orion  Population  and  how  these  properties  complement 
current  understanding  of  pre-main  sequence  stellar  evolution 
will  be  helpful  in  analyzing  the  Catalog  data  and  uncovering 
important  statistical  relationships.   The  observed  properties 
of  Orion  Population  stars  are  limited  to  the  spectra,  photom- 
etry, light  curves,  positions  and  motions.   Far  too  few  are 
known  members  of  binary  systems  to  make  the  special  tech- 
niques of  binary  studies  especially  useful  in  the  analysis 
of  the  Catalog  objects.   With  the  spectra  and  photometry  it 
is,  however,  possible  in  some  cases  to  infer  the  mass,  lumin- 
osity, composition,  dynamics  and  age  of  the  Orion  Population 
obj  ects . 

A.   OBSERVED  PROPERTIES  OF  T  TAURI  AND  T  TAURI -LIKE  STARS 

It  is  currently  believed  that  all  stars  evolve  through 
a  T  Tauri-like  stage.   Computations  of  stellar  models  show 
that  massive  objects  evolve  more  quickly  than  less  massive 
ones;  for  example,  a  star  of  one  solar  mass  will  take  about 
ten  times  as  long  to  evolve  to  the  main  sequence  as  a  star 
of  five  solar  masses.   The  main  sequence  is  a  locus  of  points 
on  a  Hertzsprung-Russell  diagram;  the  points'  positions  being 
defined  by  the  luminosities  and  effective  surface  tempera- 
tures of  stars  which  have  evolved  to  equilibrium  burning  of 
hydrogen  in  the  stellar  core.   Assuming  that  the  ratios  of 
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the  numbers  of  stars  of  different  masses,  formed  in  our 
region  of  the  galaxy,  remains  essentially  constant  with 
respect  to  time;  then  the  number  of  low  mass  pre-main 
sequence  (PMS)  stars  which  can  be  observed  is  far  more 
numerous  than  the  number  of  larger  mass  objects,  since  the 
massive  objects  spend  far  less  time  in  the  PMS  stage.   Con- 
sequently, the  largest  group  of  stars  in  the  Catalog  are 
certain  T  Tauri  stars  (39.3%)  or  presumed  so  (9.91).   The 
other  categories  with  hot  members  of  spectral  classes  B  and 
A,  late-type  low  luminosity  members  and  peculiar  stars  make 
up  only  161  of  the  Catalog.   About  3  5%  of  the  Catalog  stars 
are  not  currently  classed  at  all  but  share  some  character- 
istics.  A  thorough  discussion  of  the  T  Tauri  stars  follows, 
while  the  others  are  very  similar  and  a  brief  discussion  will 
suffice. 

1 .   Observed  Properties 

In  the  following  paragraphs  the  observed  properties 
of  T  Tauri  stars  are  discussed. 

a.   The  Spectra 

Much  more  is  now  known  about  T  Tauri  spectra 
than  was  reported  in  the  initial  establishment  of  the  class 
of  T  Tauri  stars  [Joy,  1945].   The  nature  of  the  star's 
spectrum  is  singly  the  most  important  characteristic  for 
inclusion  as  a  T  Tauri  star.   All  the  other  characteristics 
are  secondary,  although  necessary.   A  thorough  discussion 
and  delineation  of  T  Tauri  criteria  is  contained  in  a  review 
article  by  Herbig  [1962].   The  following  characteristics  arc 
observed  in  the  spectrum: 
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(1)  Emiss ion  Lines.  Numerous  emission  lines 
are  visible.  These  include  hydrogen's  (HI)  Balmer  lines, 
the  H  and  K  lines  of  once  ionized  calcium  (Ca  II),  and  in 
varying  degrees  lines  of  iron  (Fe  I,  Fe  II),  titanium 

(Ti  II)  and  forbidden  lines  of  oxygen  ([01])  and  sulphur 
([S  II]).   Only  in  late  T  Tauri  stars,  may  the  helium  (He  I) 
be  prominent.   The  hydrogen  and  calcium  lines  are  the  strong- 
est in  intensity,  followed  by  the  iron  (Fe  II)  and  titanium 
(Ti  II)  lines.1 

(2)  Underlying  Spectrum.   The  emission  lines 
of  the  previous  paragraph  are  superimposed  on  either  a  con- 
tinuous emission  spectrum  or  in  some  cases  a  nearly  normal 
photospheric  absorption  spectrum  whose  type  ranges  from  F 
to  M.   The  absorption  lines  may  be  considerably  broadened. 
In  those  stars  with  an  absorption  spectrum,  there  is  often 
a  line  of  lithium  (Li  I   6707)  which  might  legitimately  be 
included  as  a  criteria  [Kuhi,  1966].   In  general,  the  tran- 
sition from  an  absorption  spectrum  to  continuous  emission 
masking  absorption  is  most  striking  in  the  stars  with  strong 
emission  lines.   When  definite  lines  appear,  it  is  often 
possible  to  measure  doppler  shifts  separate  from  those  caused 
by  the  overall  radial  motion  of  the  star.   This  overall 
radial  motion  must  be  inferred  from  the  radial  motion  of 
associated  stars  in  the  same  cluster. 


The  Roman  Numeral  I  indicates  the  neutral  atom,  II 
indicates  the  first  ion,  etc.;  and  the  square  brackets  [  ] 
indicate  forbidden  transitions. 
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(3)  P  Cygni  Profiles.   Nearly  all  the  T  Tauri 
stars  in  which  the  absorption  lines  are  visible  have  the 
lines  arranged  so  that  there  is  an  obvious  profile  like 
that  of  P  Cygni.   This  arrangement  consists  of  an  emission 
line  symmetric  in  intensity  about  the  central  frequency  with 
an  absorption  line  bounding  the  emission  just  to  the  higher 
frequency  (blue)  side  [Ambartsumyan,  1958] .   The  exceptions 
show  an  inverse  P  Cygni  profile,  that  is,  the  absorption  is 
to  the  lower  frequency  (red)  side  of  the  emission  line 
[Walker,  1972].   These  P  Cygni  features  appear  at  the  Balmer 
lines  and  the  H  and  K  Ca  II  lines. 

(4)  Ultraviolet  Excess.   There  may  be  an  ultra- 
violet (UV)  excess  and  sometimes  a  "blue  continuum"  which 
has  been  interpreted  as  Balmer  continuum  radiation  [Strom, 
Strom,  and  Grasdalen,  1975].   This  causes  masking,  as  men- 
tioned previously,  or  as  it  is  more  commonly  called, 
"veiling."   It  may  be  characterized  by  a  steep  rise  in  the 

o 

spectral  energy  distribution  shortward  of  about  4200  A. 

b.  Polarization 

Variable  optical  polarization  as  high  as  12%  has 
been  observed.   It  is  possibly  caused  by  electron  or  dust 
scattering  [Strom,  Strom,  and  Grasdalen,  1975]. 

c.  Infrared  (IR)  Excesses 

Unusually  intense  IR  radiation  in  the  l-3y  region 
was  first  noted  by  Mendoza  [1966,  1968].   This  IR  emission 
is  clearly  greater  than  that  which  was  expected  from  fitting 
a  black  body  energy  distribution  to  the  intensities  at  visual 
wavelengths  [Rydgren,  Strom,  and  Strom,  1975]. 
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d.  Irregular  Variability 

Irregular  variability  in  brightness  was  one  of 
Joy's  original  criteria.   Continuing  analysis  in  the  past 
thirty  years  has  failed  to  uncover  any  underlying  periodicity, 
except  for  a  very  few  stars  in  which  the  time  of  independent 
characteristic  cycles  might  have  been  found..  The  range  of 
variability  extends  from  differences  of  0.2  magnitudes  to 
5.8  magnitudes.   The  mean  and  median  range  is  about  two 
magnitudes.   The  periods  range  from  several  days  to  several 
hundred  days  and  more.   There  is  variability  in  the  infrared 
also  [Mendoza,  1968],  but  which  does  not  necessarily  follow 
the  changes  in  the  visual  wavelengths  [Cohen,  1973a]. 

e.  Association  with  Nebulosity 

All  of  the  stars  appear  close  to  nebulosity. 
Of  the  T  Tauri  stars,  13.8%  illuminate  or  excite  bright 
nebulosity.   Herbig  [1962]  has  shown  that  it  is  very  unlikely 
that  these  stars  are  just  ones  which  have  been  trapped  by  the 
nebulous  clouds.   The  nebulosity  tends  to  be  heavy  and  total 
obscuration  of  other  Orion  Population  objects  is  likely 
[Strom,  Strom,  and  Grasdalen,  1975].   With  IR  surveys  and 
searches,  many  more  IR  sources  and  likely  T  Tauri  stars  are 
being  found.   It  must  be  noted  that  not  all  areas  of  great 
nebulosity  have  associated  T  Tauri  stars. 

f.  0  and  B  Associations 

T  Tauri  stars  do  appear  with  some  regularity  in 
concentrations  with  0  and  B  stars.  The  Orion  Nebula  is  the 
most  notable  example.   Because  of  their  short  lifetimes 
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the  0  and  B  stars  must  be  young,  and  the  observation  that 
the  T  Tauri  stars  appear  regularly  close  by  is  an  important 
one.   It  must  be  pointed  out,  however,  that  not  all  groups 
of  0  and  B  stars  have  the  associated  T  Tauri  stars,  the 
Pleiades   being  the  most  notable  example.   These  points  are 
discussed  more  fully  by  Kuhi  [1966] . 

g.   Non-T  Tauri  Stars  in  the  Orion  Population 
The  few  stars  which  cannot  be  classified  as 
T  Tauri  stars  usually  differ  in  their  hot  or  peculiar 
spectra.   But  they  usually  share  the  emission  character- 
istics, irregular  variability,  P  Cygni  profiles  and  IR 
excesses  of  the  T  Tauri  stars.   The  hot  members  of  the  Orion 
Population  and  other  non-T  Tauri  PMS  objects  are  discussed 
by  Strom,  Strom,  and  Grasdalen  [1975]. 

B.   INFERRED  PROPERTIES 

In  addition  to  its  importance  as  a  classification  cri- 
terion, the  spectrum  of  a  T  Tauri  star  yields  the  greatest 
information  about  those  properties  which  are  not  directly 
observable,  but  which  must  be  inferred  or  calculated  from 
the  observational  evidence.   This  paper  cannot  deal  in  depth 
with  the  calculation  or  even  estimation  of  all  the  character 
istics  and  dynamics  of  a  T  Tauri  star.   Only  the  currently 
held  conclusions  will  be  given. 

When  a  photospheric  absorption  spectrum  can  be  resolved 
or  when  an  intrinsic  color  can  be  obtained,  the  surface 
temperature  may  be  defined  to  about  two  significant  figures. 
If  the  distance  to  the  star  is  known  and  an  appropriate 
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correction  for  interstellar  (and  circumstcllar)  absorption 
can  be  applied,  then  an  absolute  magnitude  and  luminosity 
can  be  calculated  from  the  apparent  magnitude.   There  is  a 
relationship  between  mass  and  luminosity,  since  the  mass 
and  composition  uniquely  determine  all  the  physical  char- 
acteristics of  a  star.   But  T  Tauri  stars  may  not  be  static, 
equilibrium  objects.   For  those  few  stars  for  which  physical 
characteristics  can  be  estimated,  the  luminosity  is  greater 
than  what  is  expected  for  a  main  sequence  star.   One  might 
well  conclude  that  the  radius  is  greater  than  a  main  sequence 
object  of  equivalent  mass. 

These  relationships  discussed  above  can  be  described  by 
a  two-dimensional  plot  called  the  Hert zsprung-Russell  (H-R) 
diagram.   The  abscissa  of  such  a  plot  is  appropriate  units 
of  either  spectral  class,  color  or  temperature;  the  ordinate 
is  absolute  magnitude  or  luminosity.   Radius  might  be  added 
as  a  third  non-orthogonal  coordinate.   Historical  paths  of 
stars  of  a  particular  mass  can  be  plotted  on  the  H-R  diagram. 
Age  then  is  a  parameter.   In  addition,  particular  studies  of 
binary  systems  may  yield  and/or  confirm  information  about 
masses  and  radii.   Also,  there  are  indirect  ways  to  calculate 
mass  from  estimates  of  surface  gravity  when  good  spectra  are 
available. 

Good  spectra  of  T  Tauri  stars  are  exceedingly  difficult 
to  obtain  as  the  objects  are  dim  and  the  spectra  are  confused 
and  often  masked.   So,  with  this  brief  background,  the 
current  beliefs  concerning  the  T  Tauri  properties  follow. 
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1 .  Physical  Characteristics 

The  range  of  spectral  classes  has  been  discussed. 
Consequently,  the  range  of  temperatures  for  T  Tauri  stars 
follows  directly,  from  about  6600°K  to  2800°K  [Allen,  1973]. 
These  are  effective  black  body  temperatures.   The  intrinsic 
color  (B  -  V)   likewise  ranges  from  0.4  to  1.6  magnitudes. 

The  radii  are  not  at  all  well  determined,  but  aver- 
age in  the  neighborhood  of  several  solar  radii  or  greater. 

The  mass  of  the  T  Tauri  stars  probably  range  from 
0.2  to  3  solar  masses,  but  there  is  occasional  excellent 
evidence  for  larger  mass  T  Tauri  stars  [Strom,  Strom,  and 
Grasdalen,  1975].   There  are  certainly  stars  of  clearly 
greater  mass  which  exhibit  T  Tauri-like  emission  and 
variability. 

The  question  remains  as  to  how  the  mass  is  distrib- 
uted.  Because  of  the  emission  lines,  UV  continuum  and  IR 
excesses,  the  presence  of  a  circums tellar  shell  seems  likely 
This  shell  is  energized  by  a  stellar  core.   The  whole  struc- 
ture is  complicated  though,  by  rotation,  mass  loss  and  mass 
infall  discussed  below. 

2 .  Dynamics 

The  spectrum  indicates  mass  movement  and  the  rapid 
changes  in  brightness  seem  to  confirm  that  the  T  Tauri 
structure  must  be  a  dynamic  one,  though  it  is  by  no  means 
clear  whether  dynamics  of  the  envelope  or  of  the  photosphere 
are  responsible. 
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a.  Rotation 

Analysis  of  the  broadened  photospheric  absorp- 
tion lines  [Herbig,  1957]  indicates  rotation  of  the  star. 
In  comparison  with  other  late  main  sequence  stars,  the  T 
Tauri  stars  have  substantially  larger  angular  velocities 
[Rydgren,  Strom,  and  Strom,  1975].   These  velocities  range 
from  20  to  100  km/sec.   It  is  possible  that  the  broadening 
is  enhanced  by  the  expanding  circumstellar  envelope  [Strom, 
Strom,  and  Grasdalen,  1975]. 

b.  Mass  Loss  and  Infall 

The  evidence  regarding  radial  motion  of  the  T 
Tauri  stars'  mass  sometimes  seems  contradictory.   Most  of 
the  T  Tauri  stars  show  the  P  Cygni  profile.   This  indicates 
an  outflow  of  mass.   This  loss  of  mass  seems  to  be  one  in 
which  the  mass  totally  escapes  the  star.   It  has  been  esti- 
mated to  be  as  great  as  3.7  x  10    solar  masses  per  year 
[Kuhi,  1964]. 

Since  T  Tauri  stars  are  supposedly  young,  newly 
contracted  stars,  this  mass  loss  can  be  disquieting.   It  is 
believed  that  the  mass  escapes  because  there  has  not  been 
observational  evidence  of  any  of  the  ejected  matter  return- 
ing.  Prentice  [1973]  has  suggested  that  the  observed  mass 
ejection  is  a  convective  overshoot  from  the  photosphere; 
but  the  mass  is  slowed  and  the  infall  of  material  at  about 
10  km/sec  would  be  undetectable  against  the  normal  motions 
of  the  atmosphere.   It  is  believed  by  some  that  mass  loss 
is  a  major  process  in  even  very  young  stars  [Rydgren,  Strom, 
and  Strom,  1975] . 

22 


As  mentioned  previously,  some  stars  have  been 
observed  with  inverse  P  Cygni  profiles,  indicating  a  rapid 
infall  of  material  [Walker,  1972].   However,  other  inter- 
pretations regarding  the  intensity  of  various  emissions  are 
possible  [Strom,  Strom,  and  Grasdalcn,  1975].   The  weight 
of  evidence  supports  mass  ejection  for  most  T  Tauri  stars; 
and  while  some  infall  is  possible,  it  is  probably  not  part 
of  the  stellar  contraction  or  accretion  processes  as  has 
been  suggested  [Rydgren,  Strom,  and  Strom,  1975].   One  would 
expect  to  observe  spectroscopic  evidence  of  stronger  and  more 
frequent  infall  in  this  group  of  stars,  if  they  were  truly 
in  an  accretion  phase.   In  other  words,  occasional  inverse 
P  Cygni  profiles  alone  cannot  be  viewed  as  spectroscopic 
confirmation  of  an  accretion  phase. 

The  mechanics  of  mass  loss  remain  a  mystery.  It 

appears  that  the  youngest  stars  have  the  greatest  mass-loss 

rates  [Rydgren,  Strom,  and  Strom,  1975].   If  emission  line 

intensity  can  be  correlated  to  mass-loss,  then  the  observed 

fact,  that  emission  line  intensity  tends  to  become  less  for 

older  groups  of  T  Tauri  stars  [Cohen,  1974] ,  may  be  an 

important  clue  for  a  model.   Just  how  the  mass  moves  in 

connection  with  the  contraction  and  how  rotation  is  connected 

is  unknown. 

3 .   Lithium  Abundance 

Herbig  [1962]  reviews  the  work  in  the  field  of  the 

lithium  abundance  in  the  T  Tauri  stars'  composition.   The 

abundance  is  about  two  orders  of  magnitude  greater  than  the 

sun. 
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In  stars  with  hydrogen  burning  lithium  is  destroyed 
rapidly.   If  T  Tauri  stars  are  indeed  PMS  objects,  then  the 
Galactic  abundance  of  lithium  may  be  indicated.   However, 
schemes  for  the  nucleosynthesis  of  lithium  have  been  sug- 
gested.  The  actual  cause  of  this  high  abundance  is  unknown. 

C.   PRE-MAIN  SEQUENCE  (PMS)  STELLAR  EVOLUTION 

Since  it  is  now  generally  accepted  that  the  T  Tauri 
phenomenon  is  a  stage  of  stellar  evolution  through  which 
at  least  some  and  possibly  all  stars  pass,  it  is  instructive 
to  discuss  that  portion  of  a  star's  life  prior  to  the  stage 
when  it  begins  its  equilibrium  hydrogen  core  burning,  the 
main  sequence,  and  to  discuss  how  T  Tauri  stars  can  be 
assigned  to  the  PMS  set  as  support  for  a  theory  of  stellar 
evolution. 

1 .   Birth  of  a  Star 

There  are  regions  of  the  galaxy,  observationally 
close  to  earth,  in  which  molecular  hydrogen  and  dust  (only 
about  1%    of  the  interstellar  medium)  have  a  density  great 
enough  that  radiation  from  more  distant  stars  is  completely 
obscured.   These  clouds  are  pierced  by  lines  of  magnetic 
flux  and  molded  by  large  scale  gravitational  instabilities. 
They  are  bombarded  by  high  energy  cosmic  rays  and  accelerated 
by  various  Alfven  waves  in  complicated  schemes  of  coupling 
and  decoupling.   In  a  set  of  conditions,  not  well  understood, 
in  these  clouds  pockets  of  higher  density  may  be  formed.   If 
the  internal  energy  is  not  too  great,  the  sel f -gravitation 
of  the  gas  and  dust  may  be  great  enough  to  accelerate  the 
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matter  towards  a  center  of  mass.   This  is  the  start  of  a 
collapse.   By  the  Virial  Theorm  half  the  potential  energy 
released  in  the  collapse  goes  into  internal  heating  of  the 
system,  and  the  rest  must  be  radiated  from  the  system.   If 
the  mass  is  great  enough  and  the  opacity  of  the  cloud 
material  is  such  that  excess  energy  can  be  radiated  away, 
then  the  gravitational  forces  can  overcome  thermal  motions 
which  would  dissociate  the  dense  pocket,  and  further  accre- 
tion of  the  surrounding  cloud  can  continue.   Once  there  is 
a  gravitationally  bound  system,  a  condition  of  dynamic 
hydrostatic  equilibrium  exists;  and  it  is  possible  to  model 
the  system  mathematically.   This  set  of  first  order  differ- 
ential equations  relate  the  pressure,  mass,  density,  luminosity 
and  temperature  as  functions  of  a  radial  distance  from  the 
center  of  mass.   With  knowledge  or  assumption  about  the  local 
thermodynamic  state,  solutions  can  be  found.   Usually  these 
solutions  are  carried  out  in  terms  of  mass  rather  than  radius 
and  solved  numerically.   If  one  carries  out  a  series  of  such 
calculations  for  a  particular  total  stellar  mass,  then  a 
complete  history  of  the  star's  collapse,  life  and  eventual 
death  can  be  calculated.   This  history  can  be  plotted  on  the 
H-R  diagram,  showing  the  sequence  of  luminosities  versus 
effective  surface  temperatures.   Such  a  history  is  called 
a  track.   Age  is  a  parameter  along  such  a  track. 

There  are  two  distinct  tracks  currently  accepted  as 
possibilities  for  the  PMS  history.   Both  sets  of  calculations 
require  the  simplifying  assumptions  of  no  rotation,  no 
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associated  magnetic  field  and  a  spherical  fluid,  non- 
turbulent  mass.   The  earlier  set  of  calculations  led  to  what 
are  called  the  Hayashi  tracks.   These  tracks  require  the 
additional  simplifying  assumptions  of  a  roughly  homologous 
collapse  and  a  polytropic  density  distribution.   Further,  no 
provision  is  made  for  mass  loss.   More  strenuous  calculations 
yield  the  Larson  tracks.   These  different  tracks  are  discussed 
below;  however,  only  the  PMS  tracks  are  discussed,  not  the 
main  sequence  and  post-main  sequence  evolution.   Regardless 
of  model,  the  tracks  differ  greatly  depending  on  the  initial 
mass  and  the  initial  composition,  so  only  a  general  comparison 
can  be  undertaken  here,  using  a  star  of  about  a  solar  mass 
to  compare  the  two  tracks. 
2 .   The  Hayashi  Track 

A  major  result  of  Hayashi-type  calculations  is  the 
establishment  of  a  forbidden  region.   This  is  an  area  in 
the  H-R  diagram  in  which  no  stellar  object  in  hydrostatic 
equilibrium  can  exist,  because  at  any  given  mass  and  given 
luminosity,  hydrostatic  equilibrium  is  not  possible  if  the 
surface  temperature  of  the  object  is  too  cool  [Clayton,  1968]. 
In  the  Hayashi  scheme,  as  explained  by  Clayton  [1968], 
the  dense  collapsing  pocket  starts  almost  in  a  free  fall. 
Eventually  temperature  and  pressure  rise,  slowing  the  col- 
lapse to  a  quasistatic  state,  one  in  which  the  partition  of 
lost  potential  energy,  half  to  kinetic  (thermal)  energy  of 
the  particles  and  half  to  radiation  production,  becomes  a 
quite  accurate  description.   Although  the  radiation  escapes 
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easily  at  first,  eventually  there  is  increasing  opacity. 
This  is  because  opacity  is  a  state  function,  determined  by 
the  changing  density  and  temperature.   While  there  is  some 
additional  heating  because  of  this,  interior  temperatures 
cannot  rise  above  about  101*  °K  until  the  molecular  hydrogen 
has  disassociated  and  ionized.   When  the  ionization  is  com- 
plete there  can  be  further  contraction  and  heating  until  a 
state  of  dynamic  hydrostatic  equilibrium  exists,  in  the 
homologous  collapsing  sphere. 

It  is  at  this  point  that  the  Hayashi  track  begins. 
A  Hayashi  track  is  shown  in  Figure  28.   As  was  previously 
stated,  since  the  surface  temperature  cannot  be  too  cool, 
the  star's  appearance  must  take  place  at  the  near  vertical 
boundary  of  the  forbidden  region  in  the  H-R  diagram.   The 
surface  temperature  is  high  and  contraction  is  still  in 
progress,  and  the  radius  is  large.   This  must  result  in  a 
large  luminosity.   For  a  star  of  one  solar  mass,  its  starting 
luminosity  is  about  16  L~  and  effective  temperature  about 
3900°K.   The  3900°K  is  the  coordinate  of  the  boundary  of 
the  one  solar  mass  forbidden  region.   In  the  Hayashi  scheme, 
no  star  of  one  solar  mass  (with  the  same  composition)  can 
be  cooler. 

But  the  star  is  still  gravitationally  unstable  and 
collapsing,  a  process  which  continues  for  about  10   years. 
With  the  calculated  high  opacity,  the  star  must  be  fully 
convective  to  transmit  radiation  at  the  required  rate  to 
maintain  quasistatic  hydrostatic  equilibrium. 
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Subsequent  calculations  show  the  star  "moving"  nearly 
vertically  down  the  edge  of  the  forbidden  region  on  a  track 
of  only  slightly  increasing  surface  temperature  but  with  a 
dramatic  decrease  in  luminosity.   But  the  resulting  contrac- 
tion increases  central  temperature  and  changes  the  central 
opacity.   The  evolution  slows  and  a  core  in  radiative  equi- 
librium grows.   Eventually  nearly  the  entire  star  is  in 
radiative  equilibrium  and  the  main  sequence  size  is  reached, 
but  by  now  the  central  temperature  and  density  is  great 
enough  for  thermonuclear  reactions  to  occur,  supplying 
radiant  energy  at  the  same  rate  it  is  lost  at  the  surface. 
At  this  time  the  star  is  positioned  on  the  main  sequence. 
During  these  final  stages  in  which  radiative  transfer  over- 
comes the  convective  phase,  the  decreasing  luminosity  ceases 
and  the  surface  temperature  increases.   (There  is  a  small 
increase  in  luminosity  but  the  track  is  essentially  hori- 
zontal.)  Since  a  one  solar  mass  star  was  picked,  its  surface 
effective  temperature  is  about  5770°K  with  an  absolute 
bolometric  magnitude  of  +4.75. 

For  stars  of  greater  mass,  the  forbidden  region  has 
its  border  at  slightly  greater  temperatures  and  the  convec- 
tive phase  takes  a  progressively  smaller  portion  of  the 
evolutionary  time,  until  for  stars  of  the  greatest  observed 
masses  (>    15  M~),  there  is  essentially  no  convective  phase 
at  all.   In  these  cases  the  stars  move  on  a  track  of 
increasing  temperature  with  about  equal  luminosity. 
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3 .   The  Larson  Track 

The  Larson  calculation  does  not  take  the  Hayashi 
assumption  of  a  homologous  contraction  nor  the  ensuing 
polytropic  density  distribution;  but  the  assumptions  of 
spherical  symmetry,  no  rotation,  no  magnetic  fields  and 
negligible  internal  turbulent  motions  are  retained.   For 
further  detail,  the  reader  should  consult  Larson  [1969  and 
1972],  and  the  review  articles  [Strom  and  Strom,  1973]  and 
[Larson,  1973].   These  sources  were  used  for  this  summary. 

The  Larson  scheme  starts  essentially  like  the 
Hayashi  with  a  cloud  collapsing  in  free  fall,  but  not 
homologously .   The  interior  density  increases  at  the  expense 
of  the  outer  portions  whose  density  decreases.   Collapse  in 
the  outer  regions  is  "significantly  retarded  from  a  free 
fall."   This  non- homologous  collapse  results  in  a  peaked 
density  distribution,  the  greatest  density  at  the  center. 
But  furthermore,  the  collapse  becomes  most  rapid  at  the 
center  too.   The  result  is  major  changes  in  density  in  a 
decreasing  region  about  the  center  while  little  happens  in 
the  outer  regions.   Eventually,  the  central  region  becomes 
opaque  and,  with  the  rapidly  increasing  central  pressure 
and  temperature,  the  collapse  is  slowed  and  finally  stopped. 
This  core,  in  which  collapse  has  ceased,  is  in  sharp  physical 
contrast  with  an  intermediate  region  in  which  the  molecular 
hydrogen  continues  to  free  fall.   This  abrupt  change  results 
in  a  shock  front.   Infalling  material  results  in  continued 
core  collapse  subsequent  to  a  series  of  unimportant 
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rebound  oscillations.   Interestingly,  the  shock  front  which 
bounds  the  core  moves  with  the  core  surface  inward  toward 
the  center,  even  though  it  bounds  increasing  mass.   Now,  in 
this  first  core  a  second  core  is  developed  in  just  the  same 
way  except  the  molecular  hydrogen  is  dissociated  in  the 
process,  and  this  core  is  in  hydrostatic  equilibrium.   A 
second  shock  front  is  formed  but  within  the  first  core  and 
the  first  shock  front. 

A  rather  complicated  thermodynamic  condition  now 
exists  at  the  second  shock  front,  the  result  being  an  in- 
crease in  the  radius  of  the  second  core.   The  first  shock 
front  dies  away  as  the  expansion  takes  place. 

Infalling  of  matter  continues  at  the  shock  front, 
but  there  is  a  point  when  the  opacity  of  this  gas  and  dust 
drops  so  that  radiant  energy  can  be  transported  to  the  out- 
side, eventually  escaping  into  interstellar  space.   The 
result  to  an  outside  observer  would  be  a  rapid  increase  from 
almost  no  luminosity  to  about  10  times  the  luminosity  of  the 
sun.   This  is  the  start  of  the  Larson  track  for  the  luminous 
object.   A  Larson  track  is  shown  in  Figure  29.   Were  the  core 
visible,  it  would  have  a  radius  about  12  solar  radii  and  an 
effective  temperature  of  about  3200°K  for  a  single  solar  mass 
object.   But  this  is  within  the  Hayashi  forbidden  region, 
and  part  of  the  subsequent  evolution  will  carry  the  star  to 
even  lower  effective  temperatures! 

To  review  thus  far,  the  Larson  calculation  results 
in  the  formation  of  a  small  stellar  core  which  expands  while 
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accreting  additional  material;  whereas,  in  the  Hayashi  model 
there  was  a  homologous  overall  contraction  of  a  single  cloud 

Because  of  the  energy  loss,  the  core  begins  to  con- 
tract, and  a  convection  zone  starts  at  the  core's  surface, 
growing  inward.   Core  contraction  speeds  up,  fed  by  both 
radiative  and  convective  energy  losses.   The  smaller  radius 
results  in  a  lower  luminosity.   But  the  infall  of  material 
continues  and  the  further  evolution  guides  the  track  to 
greater  luminosities  and  greater  surface  temperatures,  even 
though  in  this  phase  the  radius  remains  about  constant.   In 
this  phase,  there  is  a  significant  increase  in  the  core 
mass,  and  the  track  moves  out  of  the  forbidden  region.   When 
about  half  the  total  cloud  mass  available  is  in  the  core, 
the  final  parts  of  the  track  and  movement  to  the  main  se- 
quence is  started.   At  this  point,  the  maximum  effective 
surface  temperature  has  been  reached,  about  8300°K,  and  the 
maximum  luminosity,  about  30  times  that  of  the  sun.   As  the 
remainder  of  the  cloud  mass  accretes,  the  core,  seeking  a 
resultant  static  equilibrium  state,  moves  back  to  a  lower 
temperature  and  smaller  luminosity.   As  the  outside  mass 
is  depleted,  the  kinetic  energy  of  the  infalling  material 
at  the  shock  front  is  no  longer  a  contributing  luminosity 
source.   There  is,  at  this  point,  little  infalling  material, 
but  there  is  convective  transport  in  the  outer  half  of  the 
core.   What  exists,  then,  is  a  stellar  object  at  the  edge 
of  the  forbidden  region  with  a  central  radiative  core  and  a 
convective  outer  core,  the  whole  object  in  final  collapse. 
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It  is,  in  fact,  on  the  lower  part  of  its  Hayashi  track  but 
with  a  much  smaller  radius  and  luminosity  than  was  predicted 
in  the  Hayashi  scheme.   It  did  not  reach  that  point,  however, 
in  the  same  manner  as  the  original  Hayashi  track.   From  this 
point  the  star  follows  the  Hayashi  track  to  the  main  sequence. 
In  the  case  of  protostellar  clouds  of  two  stellar  masses  and 
greater,  the  Larson  calculations  show  that  no  convective 
Hayashi  phase  is  reached  at  all. 

The  importance  of  the  Larson  scheme  lies  in  its  pre- 
diction of  a  non-homologous  collapse  resulting  in  a  stellar 
core  with  a  surrounding,  infalling  circumstellar  cloud. 
4 .   The  Evolutionary  State  of  T  Tauri  Stars 

The  question  remains  to  be  answered:  just  where  do 
T  Tauri  stars  fit  onto  an  evolutionary  track?   The  observa- 
tional clues  are:  a  stellar  object  with  visible  and  infrared 
radiation  (in  excess),  association  with  nebulosity  and  in 
some  cases  known  young  0  and  B  stars,  circumstellar  envelopes 
with  emission  and  blue  continuum,  large  rotational  velocities, 
and  mass  ejection  and/or  inflow.   Additionally,  the  luminosity 
and  temperature  can  be  found  for  a  few  representative  stars. 
These  coordinates  plot  above  the  main  sequence.   Unfortu- 
nately, their  masses  are  not  known. 

There  appear  to  be  two  main  positions  on  the  question 
of  the  T  Tauri  stars'  evolutionary  state.   The  one  delineated 
by  Rydgren,  Strom,  and  Strom  [1975]  takes  exception  to  the 
explanation  given  by  Larson  [1972]. 
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Larson  [1972]  argues  that  the  T  Tauri  stars  are 
probably  "newly  formed  stars,"  this  name  being  given  to  the 
protostar  when  it  readies  the  end  of  that  portion  of  the 
Larson  track  beginning  with  the  maximum  temperature  and 
luminosity  (half  the  mass  in  the  core)  and  ending  just  at 
the  start  of  the  Hayashi  convection  phase  (nearly  all  the 
mass  in  the  core) ,  recognizing  that  for  M  ^   2M  fl  there  is 
no  Hayashi  convection  phase.   Larson  points  out  that  the 
predicted  locus  of  "newly  formed  stars"  lies  near,  but  below, 
the  isochrone  for  106  years.   Larson's  explanation  might  fit 
some  of  the  evidence,  except  the  plotted  H-R  positions  and 
the  lack  of  observed  inverse  P  Cygni  profiles. 

This  last  fact  is  the  crux  of  Rydgren,  Strom,  and 
Strom's  [1975]  disagreement.   The  T  Tauri  stars  (a  Taurus- 
Auriga  group  and  a  p  Ophiuchus  group) ,  which  they  have 
plotted,  fall  substantially  above  the  106  year  isochrone. 
Further,  they  argue  that  optically  thick  circumstellar  dust 
shells  are  rare  among  the  T  Tauri  stars,  and  thirdly,  they 
point  out  that  the  inverse  P  Cygni  profile  is  observed  in 
only  one  of  their  sample  stars  and  that  star  has  the  lowest 
luminosity  of  the  group.   They  conclude  that  there  is  no 
"final  accretion  phase"  and  that  the  T  Tauri  stars  are 
actually  on  radiative  equilibrium  tracks  approaching  the 
main  sequence. 

While  Rydgren,  Strom,  and  Strom  may  be  correct, 
their  arguments  do  not  necessarily  refute  Larson's  thesis. 
First,  the  masses  are  not  well  known,  so  comparison  of  their 
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locus  with  Larson's  locus  drawn  for  specified  masses  is 
hardly  conclusive.   Furthermore,  the  two  groups  of  Rydgren, 
Strom,  and  Strom  have  about  the  same  locus,  yet  the  p  Oph 
group  is  known  to  be  older.   Second,  Larson's  locus  would 
change  if  rotational  effects  were  included.   This  is  a  par- 
ticularly important  point  since  the  angular  velocities  for 
T  Tauri  stars  seem  to  be  high.   The  effect  of  rotation  might 
very  well  move  the  Larson  locus  above  the  106  year  isochrome. 
Third,  there  is  no  requirement  in  the  Larson  track  for  an 
optically  thick  circumstellar  shell  at  the  time  the  "newly 
formed  stars"  appear;  rather,  if  the  star  is  appearing 
(visually)  the  shell  must  be  becoming  optically  thin.  Fourth, 
while  it  is  true  that  the  P  Cygni  profile  is  observed  with 
far  more  frequency  than  the  inverse  profile,  the  evidence 
is  not  conclusive  that  mass  is  ejected  permanently.   A  slow 
unobservable  infall  is  possible,  as  was  previously  mentioned. 
In  any  case,  at  the  point  where  Larson's  "newly  formed  stars" 
appear  there  is  no  further  requirement  for  mass  infall  since 
"practically  all  of  the  proto- stellar  mass  has  fallen  into 
the  core"  [Larson,  1969].   Furthermore,  the  mass  loss  may 
not  be  isotropic.   However,  some  explanation  for  mass 
ejection  is  still  needed. 
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III.   PRESENTATION  OF  DATA 

For  this  research  the  data  for  the  323  stars  listed  in 
the  Catalog  required  reduction  and  coding  to  make  computer 
processing  possible.   Where  qualitative  labels  were  given 
in  the  Catalog,  numerical  values  were  assigned.   Data  pro- 
cessing and  some  of  the  statistical  analyses  were  accom- 
plished using  the  Naval  Postgraduate  School  IBM  360-67 
computer  system.   The  library  programs  of  the  Statistical 
Package  for  the  Social  Sciences  (SPSS),  revision  V  [Nie, 
Bent,  Hull,  1970]  were  used  for  the  majority  of  this  work. 
Additional  regression  analyses  were  made  with  a  program 
provided  by  Mr.  Robert  Batesole  of  BDM  Corporation  that 
allowed  for  a  search  of  the  best  of  all  possible  regressions, 
and  the  regression  analysis  routines  of  Biomedical  Computer 
Programs  (BMD)  library. 

Notwithstanding  its  name,  SPSS  provided  an  excellent 
tool  for  this  project.   It  allowed  versatile  data  selection, 
recoding  and  labeling.   Further,  its  graphics  were  easy  to 
generate  and  understand. 

In  the  paragraphs  that  follow,  a  broad  look  at  the  Cata- 
log's data  will  first  be  given,  followed  by  closer  looks  at 
specific  relationships.   No  attempt  was  made  to  look  at  every 
possible  relationship  between  variables.   Some  were  clearly 
not  related  and  others  had  so  few  cases  with  complete  data 
that  analysis  was  not  significant.   Missing  values  were  not 
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allowed  to  affect  statistics.   In  addition  to  the  character- 
istics of  the  Catalog,  two  additional  variables  were  added: 
MASPRD,  which  was  calculated  from  the  difference  between 
brightest  and  dimmest  apparent  magnitudes  (m  (MAX)  -  m  (M1N)); 
and  AMPG,  the  absolute  photographic  magnitude  (M   ),  which 
was  calculated  from  distances  and  absorptions  of  prior 
research. 

One  additional  amplifying  note  needs  to  be  made.   The 
Catalog's  stars  were  divided  into  four  rough  groups  (sub- 
files) by  presumed  physical  association  as  given  in  the 
Catalog.   The  four  subfiles  used  were  defined  as  follows: 

(1)  ORION  subfile--86  stars  identified  with  the  Orion 
Nebula. 

(2)  TAUAUR  subfile-- 58  stars  identified  with  the  Taurus- 
Auriga  region. 

(3)  NGC2264  subfile--34  stars  identified  with  NGC2264. 

(4)  OTHERS  subfile--All  those  stars  not  a  member  of  the 
other  three  subfiles.   This  subfile  includes  stars 
which  are  identified  with  specific  New  General  Catalog 
(NGC)  or  Barnard  designations  even  though  they  may  be 
close  to  the  Orion  nebula  or  Taurus -Auriga  region. 
There  are  14  5  stars  in  the  OTHERS  subfile. 

Except  where  otherwise  noted,  all  stars  are  considered 
in  the  discussions  below.   Cases  with  missing  data  are 
excluded  only  in  the  considerations  of  the  particular  var- 
iable with  missing  data. 

A.   BROAD  LOOK  AT  CHARACTERISTICS 

1 .   Position  in  Galactic  Coordinates 

The  Catalog  gives  the  galactic  coordinates  expressed 
in  the  1958  system  [Allen,  1973].   The  variable  names  GALONG 
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and  GALAT  correspond  to  I         and  fa    respectively,  now  more 

2 
usually  labeled  I    and  b  . 

The  range  of  I    is  the  full  360°  while  6  has  values 
from  -34°  to  +23.3°.   Now,  referring  to  Figure  1,  it  is 
interesting  to  note  that  all  three  of  the  major  concentra- 
tions of  Catalog  stars  (ORION,  TAUAUR  and  NGC2264)  lie  in 
the  region  155°  <  b    <    215°,  a  rather  narrow  band.   Inter- 
estingly, this  band  of  high  concentration  is  bounded  on  each 
side  by  bands  of  no  Catalog  stars  at  all.   This  presumedly 
involves  the  structure  of  a  nearby  Galactic  spiral  arm. 

The  mean  value  of  b    is  -8.7°.   Excluding  the  ORION 
subfile  the  mean  is  -5.0°.   Only  after  excluding  TAUAUR 
subfile  does  b  move  within  two  degrees  of  the  Galactic  plane 
at  1.5°.   Again,  this  is  an  indicator  of  local  Galactic 
structure . 

Histograms  of  the  two  Galactic  coordinates  are  given 
in  Figures  2  and  3. 

2 .   Apparent  Magnitudes 

Photographic  apparent  magnitudes  (m   )  are  given  for 
278  stars  (86.1%  of  the  total  number  of  Catalog  stars).   No 
attempt  was  made  to  calculate  corrections  for  the  45  stars 
whose  magnitudes  were  given  at  a  visual  wavelength. 


The  Galactic  center  (  I  =    0 ,  b    =  0)  is  in  the  direction 
of  the  constellation  Sagittarius  at  right  ascension  17n  39m  3 
and  declination  -28°  54'  (1900).   In  general,  the  value  of  £ 
(GALONG)  indicates  the  direction  towards  or  away  from  the 
Galactic  center  while  b     (GALAT)  indicates  a  measure  of  dis- 
tance above  or  below  the  galactic  plane. 
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The  278  stars  with  photographic  magnitudes  have  the 


following  characteristics 


BRTMAG  DIMMAG 


mean  13  .  7  15  .  5 

std.  dev.            1.8             1.6 

max  18  ™  0  18  T  0 

mm  7.2  8.0 

Dimmest  magnitudes  were  not  given  for  118  stars 
(36.5%).   Where  only  one  magnitude  was  given,  it  was  taken 
to  be  the  brightest. 

Histograms  of  BRTMAG  and  DIMMAG  are  shown  in  Figures 
4  and  5  for  all  stars  and  thus  include  the  visual  magnitudes 
The  histograms  are  shown  for  interest  only  and  the  error 
introduced  by  the  visual  magnitudes  is  small. 
3 .   Range  of  Apparent  Magnitudes 

The  range  of  apparent  magnitudes  was  calculated  by 
taking  the  absolute  difference  between  BRTMAG  and  DIMMAG. 
There  were  205  stars  (63.5%)  for  which  such  a  range  could 
be  calculated.   The  following  statistics  were  calculated 
for  MASPRD: 

mean  2.07 

std.  dev.         1.1 
max  6.8 

min  0.2 

The  histogram  is  shown  in  Figure  6. 


Considering  only  the  certain  and  presumed  T  Tauri 
stars,  a  mean  MASPRD  of  1.89  was  found,  lower  but  not  sig- 
nificantly so.   In  calculating  MASPRD,  no  correction  was 
made  for  the  difference  in  magnitude  systems;  that  is,  the 
difference  in  visual  magnitudes  was  considered  to  be  a 
comparable  indicator  to  that  of  photographic  magnitudes. 

4 .  Light  Curves 

Only  141  stars  had  an  assigned  light  curve  classi- 
fication. Of  these  20  had  anomalous  classifications.  The 
remaining  121  stars  are  divided  among  the  classes  as  follows: 

more  frequently  bright  than  faint 30  stars  (9.31) 

most  frequently  has  the  mean  brightness .. 3 5  stars  (10.81) 

more  frequently  faint  than  bright 31  stars  (9.6%) 

no  preference  at  any  level 25  stars  {1.1%) 

5.  The  statistics  for  the  photometry  are  given  in  Table 
II.   The  histograms  of  V,  B-V  and  U-B  are  shown  in  Figures  7, 
8  and  9  respectively. 

6.  Spectra 

Spectral  types  of  Orion  Population  stars  range  from 
a  hot  B5  (T  =15,500°K)  to  a  cool  M5  (T  =2800°K).   Seventeen 
stars  have  no  classification  at  all,  14  stars  have  vague 
classifications  (e.g.,  G-K),  97  stars  appear  to  have  contin- 
uous spectra,  and  20  stars  show  clearly  continuous  spectra. 
While  the  stars  that  appear  to  have  continuous  spectra  may 
have  some  underlying  structure,  it  remains  generally  not 
resolved.   A  histogram  of  the  rough  spectral  classifications 
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appears  in  Figure  10.   This  histogram  is  for  all  the  stars 
in  the  Orion  Population. 

Considering  only  certain  T  Tauri  or  probable  T 
Tauri  stars,  the  distribution  is  like  the  one  for  all  stars. 
There  is  only  one  star  (V380  Ori)  earlier  than  F5  and  its 
spectrum  is  slightly  uncertain,  although  its  position  as 
main  sequence  dwarf  is  assumed.   As  can  be  seen  in  Figure  11, 
the  numbers  in  each  rough  spectral  class  increase  towards 
the  later  types  with  a  maximum  number  of  16  in  spectral  class 
K5.   While  only  14  stars  have  some  sort  of  M  type  spectra,  it 
has  been  pointed  out  that  this  decrease  in  numbers  is  prob- 
ably an  observational  limit  brought  about  by  the  low  visual 
magnitudes  of  these  faint  objects  [Herbig,  1962].   Stars 
earlier  than  middle  F  are  not  classed  as  T  Tauri  stars,  since 
the  specific  criteria  of  variability  and  the  usual  quality  of 
noticeable  line  broadening  (rotation)  are  not  apparent  in 
earlier  types  [Walker,  1956].   Those  early  (Be  and  Ae)  stars 
not  fitting  the  T  Tauri  criteria  are  known  as  Herbig  Emission 
stars.   They  are  similar  in  their  emission  characteristics 
and  association  with  nebulosity.   Why  the  high  rotation  and 
irregular  variability  should  start  at  middle  F  and  extend 
through  all  the  later  stars  is  not  known.   The  lesser  mass 
probably  has  some  effect  on  the  dynamics  of  stellar  evolution 
beyond  the  schemes  of  Hayashi  and  Larson. 

For  the  sake  of  completeness,  a  summary  of  the  other 
spectral  information  must  be  pointed  out.  Two  stars,  FU  Ori 
and  CPD-76°652,  have  supergiant  luminosity  classes.   Six  stars 
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have,  or  are  close  to,  giant  luminosity  classes  and  one 
star,  AN  Ori,  has  a  clear  luminosity  class  of  IV.   There 
are  41  stars  with  dwarf  luminosity  classes;  but  for  273 
stars  (84.5°o)  of  the  Orion  Population  the  luminosity  class 
cannot  be  determined  because  the  luminosity  criteria  in  the 
spectra  are  missing. 

7 .   Emission  Line  Intensity 

The  emission  line  intensity  is  a  set  of  qualitative 
categories  by  which  the  existence,  development  and  intensity 
of  a  star's  emission  relative  to  the  continuum  can  be  defined 
This  system  was  set  up  by  Herbig  [1962] .   The  emission  line 
intensity  is  not  constant,  the  lines  are  more  intense  when 
the  star  is  brighter.   The  Catalog  only  gives  one  emission 
intensity  and  does  not  indicate  at  what  point  in  the  light 
curve  the  spectrum  was  taken.   It  might  be  assumed  that  the 
Catalog  is  slightly  weighted  towards  stars  with  brighter 
emission  line  intensities  since  these  stars  are  more  likely 
to  have  been  discovered. 

The  code  assignment  increases  with  increasing  emis- 
sion intensity,  proceeding  from  no  emission  lines  to  the 
advanced  T  Tauri-type  spectrum.   A  histogram  of  the  emission 
intensities  is  shown  in  Figure  12. 

Cohen  [ 1974]  has  described  how  younger  clusters  tend 
to  have  distributions  such  that  there  is  a  greater  number 
proportionally  of  T  Tauri  stars  having  greater  emission 
intensities.   In  older  clusters  the  number  densities  tend 
to  favor  lesser  categories  of  emission  intensities. 
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Cohen's  results  were  checked  using  the  much  larger  sample 

size  of  this  study  to  plot  the  histograms  and  calculate 

chi-squared  statistics  for  testing  his  hypotheses.   Cohen's 

results  are  supported. 

8 .   Types  of  Stars  in  the  Catalog 

Given  below  is  the  number  of  stars  assigned  to  each 

rough  specification  of  the  star's  type: 

no  information 114 

certain  T  Tauri 127 

presumed  T  Tauri 32 

hot  member  (usually  of  spectral  class  B  or  A) .  .  .  .  21 

a  late  type  star  of  low  to  moderate  luminosity  ...  20 

peculiar  object   9 

B.   ABSOLUTE  MAGNITUDES 

On  the  basis  of  thorough  literature  research  for  dis- 
tances, reddening  and  absorption,  absolute  photographic, 
magnitudes  were  calculated  for  228  stars.   The  distance 
moduli,  (m    -  M   ),  and  absolute  photographic  magnitudes 
(M   )  are  listed  in  Table  III  along  with  the  source. 
1 .   Calculation  of  Distance  Modulus 

In  general,  a  distance  and  monochromatic  extinction 
in  magnitudes  is  required  to  calculate  a  distance  modulus, 
(m  -  M) ,  from  the  well  known  relationship 

(m  -  M)  =  5  log(d)  -  5  +  A 
where:   m  is  the  apparent  magnitude: 
m  =  constant  -  2.5  log  F 

=  constant  -  2.5  log  (^-jir) 

=  constant  +  5  log  d  -  2.5  log  L   ; 
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and  where:   d  is  the  distance  from  earth  to  the  star, 
usually  in  parsecs  (pc) ,  F  is  the  monochromatic  flux,  i 
is  the  monochromatic  luminosity,  and  A  is  an  absorption 
in  magnitudes ; 

and  where:   M  is  the  absolute  magnitude  which  is  defined 
as  the  apparent  magnitude  of  a  star  at  a  standard  distance 
of  10  pc,  that  is 

M  =  constant  +  5-2.5  log  I  .     . 
2 .   Calculation  of  Absorption 

The  value  of  the  absorption  is  rather  more  compli- 
cated to  calculate  because  a  star's  light  is  not  only  dimmed 
but  reddened.   Because  longer  wavelengths  are  not  absorbed 
to  the  same  extent  as  shorter  wavelengths  it  is  possible  to 
calculate  a  ratio  (R)  of  total  to  selective  absorption  when 
photometry  and  an  intrinsic  (unreddened)  color  can  be  found. 
It  is  unfortunate  that  this  ratio  is  not  everywhere  constant 
Newly  calculated  values  were  used  wherever  available;  but 
when  R  was  unknown  the  standard  value  R  -  3  was  used.   With 
a  value  of  R,  the  absorption  A  was  calculated  from 

A  =  R  x  E, 

where  E  is  the  color  excess,  that  is,  the  difference  between 
the  observed  and  intrinsic  colors.   Colors  are  defined  as 
the  difference  in  magnitude  between  apparent  magnitudes, 
which  are  observed  at  two  different  effectiA^e  wavelengths, 
for  example ,  B  -  V  . 
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For  those  stars  for  which  no  value  of  R  or  E  could 
be  found,  a  standard  absorption  was  assumed  depending  on 
the  stars'  galactic  latitude. 

3 .  Calculation  from  Luminosity 

In  a  few  cases  absolute  magnitudes  were  determined 
directly  from  the  total  luminosity  using  the  below  rela- 
tionship and  bolometric  corrections  from  Allen  [1973]. 

M.  ,  =  4.75  -  2.5  log!  (in  solar  units), 
bol 

4 .  Results  of  Calculations 

The  formulae  above  were  used  whenever  the  data  per- 
mitted.  Special  assumptions  regarding  the  calculations  are 
contained  in  the  Notes  to  Table  III.   Except  where  noted, 
all  magnitudes  were  corrected  to  photographic  magnitudes. 
A  histogram  of  the  absolute  magnitudes  is  shown  in  Figure 

13.   The  distribution  of  stars  in  M   by  subfile  groups  is 

Pg   '         fa    i 

discussed  in  a  later  section  on  the  luminosity  function. 
The  following  statistics  were  calculated  for  M 


Pg 


mean  4  .  8 


std.  dev.         2  1  1 


max  e . a 

min  -4.1 


C.   RELATIONSHIPS  BETWEEN  CHARACTERISTICS 

There  are  many  cross  relationships  that  can  be  formed 
between  various  categories  given  in  the  Catalog.   Many  of 
these  were  studied  using  Cartesian  scatterplots  and  two-way 
frequency  tables.   Where  indicated,  second  order  effects 
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were  examined.   Individual  H-R  diagrams,  using  apparent 
magnitudes,  for  Taurus-Auriga,  Orion  and  NGC2264  were 
examined;  but  differential  reddening  and  extinction  badly- 
scattered  the  results,  and  so  the  diagrams  were  not  helpful. 
This  problem  may  still  be  a  factor  even  in  H-R  plots  using 
M   since  for  ORION  and  NGC2264  an  average  extinction  was 
used  throughout,  a  clear  probable  error  of  up  to  a  magni- 
tude, at  least  in  the  ORION  subfile.   For  "the  few  stars 
with  a  characteristic  period,  no  relationship  was  found 
between  the  characteristic  period  and  the  star's  absolute 
magnitude . 

Three  sets  of  relationships  of  particular  interest  are 
discussed  below. 

1 .   Correlation  with  the  Range  of  Apparent  Magnitudes 
with  Other  Characteristics 

When  this  research  began,  it  was  assumed  that  the 
range  in  apparent  magnitude  (MASPRD)  was  related  to  the 
dynamics  of  the  star  in  some  way;  and  that  being  so,  a 
correlation  would  be  found  between  MASPRD  and  some  other 
observable  characteristic.   In  a  thorough  analysis,  riuo 
significant  correlation  between  MASPRD  and  any  other  variable 
or  set  of  variables  could  be  found. 

MASPRD  was  plotted  on  two-dimensional  scatterplots 
against  each  of  the  following  variables  without  significant 
correlation:   GALONG,  GALAT,  BRTMAG ,  DIMMAG,  V,  BMNSV,  UMNSB, 
TEMP,  EMINT,  and  AMPG.   The  highest  correlation  was  with 
BRTMAG  (r  =  0.4).   A  two-way  frequency  distribution  was 
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prepared  for  MASPRD  and  LICRV.   There  was  no  essential 
difference  in  the  distribution  of  MASPRD  for  the  four  main 
light  curve  classes.   For  the  variable  TEMP  only  those  stars 
with  well  defined  non-continuous  spectra  were  considered, 
since  stars  with  continuous  spectra  were  given  identical 
codes.   Thus  there  was  no  correlation  between  MASPRD  and 
the  effective  photospheric  temperature. 

There  were  61  stars  that  had  values  of  BRTMAG , 
DIMMAG,  V,  BMNSV,  TEMP,  EMINT  for  all  cases.   A  program  was 
prepared  to  find  the  best  of  all  possible  regressions  of 
these  variables  with  MASPRD  as  the  dependent  variable.   A 
linear  combination  of  BRTMAG  and  V  provided  the  best  corre- 
lation; however,  the  calculation  of  an  F  statistic  confirmed 
the  poorness  of  fit.   A  stepwise  regression  was  also  run 
with  similar  results.   A  plot  of  the  residuals  versus  the 
variables  showed  no  underlying  structure.   Factor  analysis 
was  not  undertaken,  but  may  show  some  significant  tie  undis- 
covered in  the  regression  analysis.   This  possibility  is 
considered  unlikely.   At  this  time,  there  is  no  clear  rela- 
tionship between  the  spread  in  apparent  magnitudes  and  any 
other  observable  characteristic  or  set  thereof. 

This  is  particularly  disturbing  in  the  case  of  the 
lack  of  any  relationship  between  MASPRD  and  the  emission 
intensity  (EMINT),  a  scatterplot  of  which  is  shown  in  Fig. 
14.   There  are  many  indications  that  the  emission  intensity 
and  variability  are  related.   As  stated  previously,  the 
emission  lines  are  stronger  when  the  star  is  at  maximum  light 
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Some  strong  emission  line  stars  show  small  light  variations 
on  a  short  time  scale  in  addition  to  the  large  variations 
characterized  by  MASPRD,  while  weak  emission  line  stars  do 
not  show  the  small  variations  [Rydgren,  Strom,  and  Strom, 
1975].   It  has  been  asserted  by  Rydgren,  Strom,  and  Strom 
[1975]  that  "the  principal  source  of  variability  in  T  Tauri 
stars  is  the  changing  strength  of  the  envelope  emission." 
Lastly,  those  characteristics  which  delineate  PMS  objects 
should  lessen  in  degree  and  disappear  as  the  star  evolves 
to  the  main  sequence. 

Since  main  sequence  stars  have  neither  irregular 
variability  nor  emission  lines,  and  the  PMS  T  Tauri  stars 
do,  presumedly  both  emission  strength  and  variability  de- 
crease in,  at  least,  a  qualitative  way  with  increasing 
stellar  age.   But  Figure  14  shows  no  tendency  for  MASPRD 
to  change  with  EMINT.   Statistical  tests  confirm  the  visual 
observation  that  there  is  no  correlation.   No  current 
theories  explain  either  the  cause  of  the  irregular  varia- 
bility, nor  why  MASPRD  should  be  independent  of  the  other 
observed  characteristics. 

The  question  of  the  correlation  between  the  magni- 
tude of  the  variability  and  the  emission  intensity  may  be 
important  to  determining  the  true  nature  of  certain  A  and 
F  stars.   Strom,  Strom,  and  Grasdalen  [1975]  have  suggested 
that  if  the  variability  and  envelope  emission  strength  are 
correlated,  then  certain  weak  Ha  stars  may  be  small  amplitude 
variables.   No  correlation  can  be  suggested  on  the  basis  of 

this  research. 
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2  .   Hert zsprung-Russcll  Diagrams 

The  poor  results  obtained  from  plotting  apparent 

magnitudes  versus  either  the  spectral  class  or  the  color 

B-  V  have  already  been  discussed.   The  use  of  the  absolute 

magnitude  improves  these  plots  considerably,  particularly 

in  the  H-R  diagram  of  the  Taurus-Auriga  region.   The  list 

below  shows  some  of  the  H-R  diagrams  plotted  in  this 

research : 

Figure  15  -  M   vs.  TEMP  -  all  stars 
Pg 

Figure  16  -  M    vs.  TEMP  -  subfile  OTHERS 
h  Pg 

Figure  17  -  M   vs.  TEMP  -  subfile  TAUAUR 
fe  Pg 

Figure  18  -  M    vs.  TEMP  -  subfile  NGC2264 
b  Pg 

Figure  19  -  M    vs.  TEMP  -  subfile  ORION 
s  Pg 

Figure  20  -  M   vs.  TEMP  -  all  certain  or  probable 
P^  T  Tauri  stars 

Figure  21  -  M   vs.  B  -  V  -  all  stars 
b  Pg 

Figure  22  -  M   vs.  B-V  -  subfile  OTHERS 
6  Pg 

Figure  25  -  M    vs .  B  -  V  -  subfile  TAUAUR 


& 


pg 


Figure  24  -  M    vs.  B-V  -  subfile  NGC2264 
b  Pg 

Figure  25  -  M    vs.  B-V  -  subfile  ORION 
h  Pg 

Figure  26  -  M   vs.  B-V  -  all  certain  or  probable 
P£  T  Tauri  stars 

The  dwarf  main  sequence  is  drawn  on  each. 

Two  items  regarding  these  diagrams  are  causes  of 

error.   First,  most  of  the  stars  (ORION  and  NGC2264)  had 

their  absolute  magnitudes  calculated  on  the  basis  of  an 

average  extinction.   This  could  cause  an  estimated  standard 

error  of  as  much  as  1  magnitude.   However,  since  the  error 
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is  applied  to  stars  equally  in  both  directions,  the  overall 
effect  probably  does  not  bias  the  considerations  of  large 
scale  characteristics  from  the  H-R  diagrams.   On  the  other 
hand,  the  absolute  magnitudes  for  the  Taurus -Auriga  stars 
have  quite  small  standard  errors,  probably  no  greater  than 
±0.3  magnitudes.   This  is  because  much  closer  estimates  of 
individual  extinction  could  be  applied,  given  the  excellent 
photometry  and  calculations  of  Rydgren,  Strom,  and  Strom 
[1975J.   Second,  but  for  equivalent  reasons,  the  H-R  dia- 
grams with  B-V  as  the  abscissa  are  in  even  greater  error, 

because  not  only  are  the  M   ' s  in  error  as  above;  but  the 

7  Pg 

reddening  is  variable  also.   Consequently,  a  star's  hori- 
zontal position  is  not  only  offset  to  a  greater  B-V  because 
of  overall  reddening,  but  is  uncertain  by  an  amount  not 
determined.   Again,  for  recognition  of  gross  characteristics 
these  errors  are  probably  inconsequential. 

After  examination  of  the  diagrams  there  can  be  little 
doubt  that  Orion  Population  stars  usually  have  luminosities 
greater  than  main  sequence  stars  of  the  same  spectral  class, 
at  least  for  mid-F  stars  and  later.   There  are  seven  stars 
earlier  than  F5  with  peculiar  luminosities.   Particularly 
strange  are  LkHa  198  (HRC  3)  and  LHa  25  (HRC  219).   These 
stars  are  too  dim  by  about  three  magnitudes.   LHa  25  is  in 
NGC2264,  a  cluster  with  a  well  determined  distance  modulus 
[Walker,  1956].   Its  spectrum,  B8pe  +  shell,  indicates  a 
source  of  circumstellar  extinction.   The  most  probable  cause 
of  the  error  is  a  large  difference  in  the  extinction  from 
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the  average  value  assumed  for  NGC2264.   The  distance  modulus 
for  the  small  anonymous  dark  nebula  close  to  Lkikx  198  has 
been  determined  spcctroscopically  from  two  main  sequence  B 
stars  [Herbig,  I960].   LkHa  198*s  spectrum  is  an  uncertain 
A  and  so  the  necessary  absorption  to  take  the  star  to  the 
main  sequence  ranges  from  two  to  four  magnitudes.   Extinc- 
tion is  probably  the  cause  of  the  anomaly,  but  too  little 
is  known  about  this  nebula  to  draw  any  firm  conclusions. 
The  other  five  stars  may  truly  be  large  mass  PMS  objects  or 
the  absorption  has  been  badly  overestimated. 

For  the  F5  and  later  stars  the  H-R  diagrams  combined 
with  the  various  spectral  peculiarities  would  seem  to  con- 
firm their  PMS  nature.   On  the  basis  of  these  H-R  diagrams 
alone,  particular  "tracks"  cannot  be  chosen.   But  from 
Figure  15  one  can  see  that  the  cross  section  of  Catalog 
stars  covers  stages  of  evolution  from  early  visual  appearance 
to  positioning  on  the  main  sequence,  if  Larson  tracks  are 
assumed. 

3.   The  Relationship  Between  the  Light  Curve 
Classification  and  the  Emission  Intensity 

It  has  been  previously  shown  that  there  is  no  cor- 
relation between  EMINT  and  MASPRD  and  further,  that  the 
distribution  of  MASPRD  does  not  change  between  the  various 
light  curve  classifications.   It  seems  reasonable  that  if 
EMINT  and  MASPRD  are  independent  and  if  the  distribution  of 
MASPRD  for  each  LICRV  class  is  the  same  (and  possibly  random), 
that  there  might  be  some  information  gained  by  considering  a 
two-way  frequency  table  of  LICRV  and  EMINT.   This  cross 
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tabulation  is  shown  in  Table  IV.   In  making  this  table, 
intermediate  categories  of  EMINT,  that  is,  EMINT  =  1.5,  2.5, 
3.5  or  4.5,  were  counted  as  the  next  highest  value.   For 
example,  cases  with  EMINT  =  1.5  are  counted  as  cases  with 
EMINT  =2.0.   A  chi-squared  statistic  was  calculated  to  test 
the  hypotheis  that  the  distributions  of  EMINT  between  light 
curve  classifications  were  equivalent.   The  distributions 
of  EMINT  differ  at  a  991  confidence  level.   This  may  clearly 
be  seen  by  examination  of  Table  IV,  row  by  row.   The  inter- 
esting effect  is  that  those  stars  which  are  more  frequently 
bright  than  faint  have  clearly  the  weakest  emission  intensity 
Recall  that  for  individual  stars,  emission  line  strengths 
are  greatest  when  the  stars  are  brightest.   Furthermore, 
this  difference  is  still  distinct  if  a  statistical  collapse 
is  made  with  intermediate  categories  moved  to  the  next  lower 
value  or  even  split  between  values.   The  light  curve  clas- 
sification with  the  strongest  emission  intensity  is  the 
group  which  most  frequently  has  the  mean  brightness.   No 
simple  explanation  is  apparent. 

D.   THE  LUMINOSITY  FUNCTION 

The  luminosity  function  shows  the  relative  numbers  of 

stars  in  successive  intervals  of  absolute  magnitude  for  a 

given  volume  of  space.   In  Figure  27  the  histograms  of  M 

for  each  of  the  three  specific  subfiles  and  for  the 

luminosity  function  [Allen,  1973]  are  shown.   The  luminosity 

function  (M   )  gives  the  number  of  stars  per  unit  volume 
Pg  & 

within  the  magnitude  range  M-  1/2  to  M  +  1/2.   The  sharp 
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divergence  between  the  observed  results  and  the  luminosity 

function  at  about  M    =  +6  is  probably  due  in  part  to  ob- 

pg 

servational  selection.   The  faint-end  cut-off,  however, 
appears  to  be  too  sharp  to  be  explained  by  observational 
limitations  alone;  thus  the  luminosity  function  for  T  Tauri 
stars  appears  to  be  significantly  different  than  the  lumin- 
osity function  for  the  solar  neighborhood.   This  difference 
is  not  predicted  by  theory. 
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IV.   CONCLUSIONS 

The  major  results  of  this  research  were  estimates  of 
absolute  magnitude  for  228  Orion  Population  stars,  the 
calculation  of  a  variety  of  statistics  for  the  Catalog's 
categories,  the  production  of  frequency  distributions,  and 
the  creation  and  statistical  testing  of  many  cross  relation- 
ships between  the  Catalog's  variables.   This  included 
Hertzsprung-Russell  diagrams.   These  H-R  diagrams  are  quite 
accurate  particularly  in  the  Taurus -Auriga  region.   There 
is  no  Orion  Population  star  F5  or  later  with  an  anomalous 
position. 

The  Orion  Population  stars,  of  which  T  Tauri  stars  are  a 
major  subset,  are  largely  pre- main  sequence  objects  evolving 
towards  the  main  sequence.   All  current  theories  of  stellar 
evolution  predict  the  appearance  of  stars  with  H-R  diagram 
positions  that  the  T  Tauri  stars  have  been  observed  to  fill. 
The  majority  of  the  plotted  stars  fall  well  below  the  giant 
star  branch.   However,  on  the  basis  of  large  scale  consider- 
ations of  the  Catalog  data,  no  conclusive  argument  can  be 
made  for  any  particular  track,  convective  or  radiative. 

One  could  certainly  reach  this  conclusion  only  on  the 
basis  of  a  lack  of  understanding  of  the  hydrodynamics.   The 
theories  of  Larson  and  Hayashi  do  not  explain  the  mass  loss, 
variability,  or  our  results.   On  the  basis  of  this  work 
there  does  not  seem  to  be  any  connection  between  the  large 
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scale  variability  as  represented  by  MASPRD  and  any  other 
observable  characteristic.    Since  one  expects  the  directly 
observable  characteristics  to  be  linked  through  the  equations 
of  hydrodynamic  equilibrium,  clearly  the  full  dynamical  sit- 
uation cannot  yet  be  understood.   There  does  seem  to  be  a 
relationship  between  the  emission  intensity  .and  light  curve 
classification,  but  this  is  also  not  yet  explainable. 

The  sharp  cut-off  observed  in  the  luminosity  function 
indicates  a  population  for  T  Tauri  stars  differing  from  the 
solar  neighborhood.   The  cut-off  appears  in  each  of  the  major 
groups  of  T  Tauri  stars.   It  may  be  indicative  of  T  Tauri 
structure  as  yet  unexplained  in  theory. 
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TABLE  I 


CATEGORIES  OE  THE  HERBIG-RAO  CATALOG 


10 
11 
12 
13 
14 
15 
16 
17 
18 
19 
20 

21 

22 


23 


24 


Category 

HRC  Number 

Indication  of  associated  bright  nebulosity 

Star  name  or  other  designation 

Right  ascension  and  declination  (1950) 

Galactic  coordinates  (1958) 


"I" 
"fa" 


Apparent  magnitudes  (m ) : 

Brightest  magnitude 
Dimmest  magnitude 

Light  curve  classification 

Comments  regarding  the  light  curve 

Logarithm  of  a  characteristic  time 
interval  in  days,  over  which  the  star 
has  been  observed  to  undergo  a  major 
change  in  light  and  to  recover 

Mean  visual  photometric  magnitude  (V) 

Mean  value  of  the  color  B-V 

Mean  value  of  the  color  U-B 

Ultraviolet  excess  information 

Uncertainty  of  photometric  data 

Comment  regarding  variations  in  V 

Spectral  class 

Luminosity  class 

Information  of  emission  lines 

Rough  specification  of  the  type  of  star 

A  number  describing  the  general  intensity 
of  the  emission  spectrum 

Code  indicating  if  Herbig  examined 
spectrum  plate 

Indication  of  abnormally  strong 
ultraviolet  continuum 

Categories  not  in  the  Catalog 
but  Used  in  this  Project 

The  difference  between  the  brightest 
apparent  magnitude  and  the  dimmest 
apparent  magnitude,  i.e.,  DIMMAG  -  BRTMAG 

The  absolute  photographic  magnitude,  M 


Fortran  Name 

HRC 

NEB 

none 

none 

G A LONG 
GALAT 

BRTMAG 
DIMMAG 

LICRV 

LI  COM 

CRVLOG 


V 

BMNSV 

UMNSB 

UVEX 

MUNC 

IVARV 

TEMP 

L 

IE 

ITYPE 

EMINT 

H 

IUCON 


Pg 


MASPRD 


AMPG 
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TABLli  II.   PHOTOMETRY 

Photometry  of  Orion  Population 

V        B-V        U-B 


Valid  observations  137  134  120 

Missing  observations  186  189  203 

Mean  12  *?  93  1  ™  09  0  T  24 

Median  13  ?  26  1  ?  12  0  ?  23 

Std.  dev.  1.9  0.4  0.6 

Max.  16  7  94  2  ?  14  2  ?  33 

Min.  7  ?  01  0  7  14  -1  I1  32 


Photometry  of  Certain  or  Presumed  T  Tauri  Stars 

V  B^V  U-B 

Valid  observations          72  72  65 

Missing  observations        87  87  94 

Mean                    13  ?  24  1  7  14  0  7  11 

Median                  13  ?  48  1  7  16  0  ™  05 

Std.  dev.                  1.6  0.3  0.6 

Max.                     16  ?  94  1  7  80  2  ?  33 

Min.                     9  ?  24  0  ?  50  -1  ?  32 
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TABLE  III 
STAR 

.   ABSOLUTE 

PHOTOGRAPHIC  MAGNITUL' 

lES 

HRC 

LOCATION 

fm -Ml 
JVR 

M 
Pg 

REFERENCE 

NOTE 

1 

MacC  HI 2 

CEP  IV 

12.8 

4.4 

1, 

2 

2 

Lk  Ha  19  7 

ADN 

9.9 

7.1 

3 

1 

3 

Lk  Ha  198 

ADN 

9.9 

5.1 

3 

4 

MacC  H10 

CEP  IV 

12.8 

3.6 

1, 

2 

5 

MacC  H9 

CEP  IV 

12.8 

3.3 

J, 

2 

6 

Lk  Ha  200 

NGC  225 

10.2 

3.8 

1, 

2 

7 

Lk  Ha  201 

NGC  225 

10.2 

3.8 

1, 

2 

17 

L  Ha  92 

IC  348, B4 

9.2 

7.3 

4 

18 

L  Ha  97 

IC  348, B4 

9.2 

8.8 

4 

20 

Lk  Ha  330 

IC  348 

9.2 

3.8 

4 

40 

Lk  Ha  101 

NGC  1579 

13.9 

3.1 

5 

84 

CO  Ori 

B  2  25 

8.1 

3.  5 

6 

85 

GIV  Ori 

B  225 

8.1 

2.6 

6 

86 

V649  Ori 

B  225 

8.1 

4.9 

6 

164 

V380  Ori 

NGC  1999 

10.  5 

0.3 

7 

183 

V625  Ori 

B  35 

12.  6 

2.2 

7 

2 

184 

V630  Ori 

B  35 

12.6 

2.5 

7 

2 

185 

V631  Ori 

B  35 

12.6 

2.4 

7 

2 

186 

FU  Ori 

B  35 

12.6 

-2.4 

7 

187 

Lk  Ha  314 

NGC  2068 

8.6 

6.4 

8 

3 

207 

R  Mon 

NGC  2  2  61 



-2.  7 

7 

251 

RU  Lup 

ADN 



4.0 

9 

4 

255 

Do-Ar  9 

B  42 

9.6 

3.4 

10 

5 

256 

Haro  1-1 

Near  B  42 

8.9 

5.  5 

10 
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TABLE  III  (continued) 


HRC 

STAR 

LOCATION 

(m-M) 

Pg 

M 
Pg 

REFERENCE 

NOTE. 

257 

Haro  1-4 

B 

42 

10.3 

3.8 

10 

258 

V8  5  2  Oph 

B 

42 

9.6 

5.2 

10 

5 

259 

S-R  4 

B 

42 

9.7 

3.3 

10 

260 

Do-Ar  22 

B 

42 

9.6 

3.4 

10 

5 

261 

Haro  1-8 

B 

42 

10.0 

5.5 

10 

262 

S-R  24 

B 

42 

9.6 

5.4 

10 

5 

263 

S-R  12 

B 

42 

9.6 

3.9 

10 

5 

264 

S-R  9 

B 

42 

8.9 

3.6 

10 

2  65 

S-R  10 

B 

42 

8.7 

6.3 

10 

266 

V853  Oph 

B 

42 

,44 

9.6 

3.  7 

10 

5 

267 

Haro  1-14 

B 

42 

9.6 

6.2 

10 

5 

268 

Haro  1-16 

B 

44 

10.6 

5.3 

10 

277 

Lk  Ha  122 

NGC 

6514 

12.2 

3.8 

4 

278 

V1752  Sgr 

NGC 

6523 

12.2 

4.7 

2 

279 

Lk  Ha  102 

NGC 

6523 

12.2 

2.8 

2 

280 

SV  Sgr 

NGC 

6523 

12.6 

2.1 

2 

281 

Lk  Ha  118 

NGC 

6523 

12.2 

-0.2 

2 

288 

R  CrA 

NGC 

6729 

13.1 

-1.9 

7 

290 

T  CrA 

NGC 

6729 

10.2 

2.2 

7 

299 

V521  Cyg 

NGC 

7000 

14.7 

0.8 

7 

6 

300 

V1057  Cyg 

NGC 

7000 

14.  7 

-4.1 

7 

301 

Lk  Ha  191 

NGC 

7000 

14.7 

-0.7 

7 

6 

308 

Lk  Ha  349 

IC  1396 

12.7 

2.3 

4 

310 

BD+46°3471 

IC  5146 

11.7 

-1.3 

1 

311 

Lk  Ha  245 

IC  5146 

11.7 

4.9 

1 
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TABLE  III  (continued) 

HRC     STAR      LOCATION     (m-M)      M     REFERENCE    NOTE 
Pg     Pg 

312  Lk  Ha  257 

314  Lk  Ha  350 

319  MacC  H3      Cep  IV       12.8      2.7     1  7 

320  MacC  113     Cep  IV       12.8     3.1     1         7 

321  Lk  Ha  259    Cep  IV       12.8      3.7     1  7 

322  MacC  H5      Cep  IV       12.8      4.6     1  7 

323  MacC  H18     Cep  IV       12.8      2.6     1  7 

23  FM  Tau      Tau-Aur,       9.9     4.8    10         8 

24  FN  Tau 

2  5  CW  Tau 

26  FP  Tau 

27  CX  Tau 

28  CY  Tau       Tau-Aur,       9.7      5.0    10  9 

B27 

29  V410  Tau     Tau-Aur,       9.7      2.3    10  9 

30  DD  Tau 

31  CZ  Tau 

32  BP  Tau 

33  DE  Tau 

3  4  RY  Tau 

3  5  T  Tau 

36  DF  Tau      Tau-Aur        9.0      3.5    10  10 

37  DG  Tau      Tau-Aur,       9.0      3.5    10  10 


IC  5146 

11.7 

2.1 

1 

Cep  III 

12.  8 

3.2 

1 

Cep  IV 

12.8 

2.7 

1 

Cep  IV 

12.8 

3.1 

1 

Cep  IV 

12.8 

3.7 

1 

Cep  IV 

12.8 

4.6 

1 

Cep  IV 

12.8 

2.6 

1 

Tau-Aur, 
B209 

9.9 

4.8 

10 

Tau-Aur, 
B209 

9.7 

5.9 

10 

Tau-Aur, 
B209 

11.6 

4.4 

10 

Tau-Aur 

6.4 

8.6 

10 

Tau-Aur 

6.  5 

8.6 

10 

Tau-Aur, 
B27 

9.7 

5.0 

10 

Tau-Aur, 
B7 

9.7 

2.3 

10 

Tau-Aur, 
BIO,  B7 

10.2 

5.1 

10 

Tau-Aur, 
B10,  B7 

9.1 

7.6 

10 

Tau-Aur, 
near  B7 

7.9 

4.9 

10 

Tau-Aur 

8.8 

5.  5 

10 

Tau-Aur, 
B214 

8.4 

3.3 

10 

Tau-Aur 

9.0 

2.5 

10 

Tau-Aur 

9.0 

3.  5 

10 

Tau-Aur, 
B217 

9.0 

3.5 

10 
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TABLE  III  (continued) 

HRC     STAR       LOCATION     (m-M)      M     REFERENCE    NOTE 

Pg     Pg 


38 

DH 

Tau 

Tau-Aur , 
B19 

7.3 

7.7 

10 

11 

39 

DI 

Tau 

Tau-Aur , 
B19 

7.3 

7.0 

10 

41 

IQ 

Tau 

Tau-Aur, 
B218 

9.0 

6.3 

10 

10 

43 

UX 

Tau 

A 

Tau-Aur 

8.5 

3.7 

10 

44 

FX 

Tau 

Tau-Aur , 
B18 

8.7 

6.1 

10 

12 

45 

DK 

Tau 

Tau-Aur , 
near  B19 

9.75 

3.7 

10 

46 

zz 

Tau 

Tau-Aur , 
B18 

8.7 

5.5 

10 

12 

47 

Lk 

Ha  331 

Tau-Aur , 
B18 

8.7 

7.8 

10 

12 

48 

HK 

Tau 

Tau-Aur, 
B18 

8.7 

7.6 

10 

12 

49 

HL 

Tau 

Tau-Aur 

13.8 

1.2 

10 

12 

50 

XZ 

Tau 

Tau-Aur 

13.8 

1.7 

10 

13 

51 

Lk 

Ha  266 

Tau-Aur 

10.0 

5.0 

10 

53 

UZ 

Tau 

P 

Tau-Aur, 
B19 

9.0 

5.0 

10 

10 

54 

GG 

Tau 

Tau-Aur 

9.0 

4.7 

10 

55 

GH 

Tau 

Tau-Aur, 
B18 

8.7 

6.0 

10 

12 

56 

GI 

Tau 

Tau-Aur , 
B18 

9.4 

5.7 

10 

57 

GK 

Tau 

Tau-Aur , 
B18 

9.4 

4.7 

10 

58 

DL 

Tau 

Tau-Aur, 
near  B18 

10.0 

4.  5 

10 

59 

IS 

Tau 

Tau-Aur, 
B19 

9.0 

7.  5 

10 

10 

60 

HN 

Tau 

Tau-Aur 

10.8 

3.7 

10 

61 

CI 

Tau 

Tau-Aur 

10.2 

4.0 

10 

13 

60 


TABLE    III    (continued) 


HRC 

STAR 

LOCATION 

(m-M) 

M 
Pg 

REFERENCE 

NOTE 

62 

DM 

Tau 

Tau-Aur 

6.9 

7.8 

10 

63 

AA 

Tau 

Tau-Aur, 
B18 

8.4 

5.2 

10 

64 

HO 

Tau 

Tau-Aur 

8.8 

7.4 

10 

13 

65 

DN 

Tau 

Tau-Aur , 
B18 

7.5 

6.2 

10 

66 

HP 

Tau 

Tau-Aur 

8.8 

7.2 

10 

13 

67 

DO 

Tau 

Tau-Aur , 
B22 

11.7 

3.4 

10 

68 

VY 

Tau 

Tau-Aur 

7.5 

5.8 

10 

13 

69 

Lk 

Ha  332 

Tau-Aur , 
B22 

10.  2 

6.3 

10 

13 

70 

DP 

Tau 

Tau-Aur, 
B22 

10.2 

5.0 

10 

13 

71 

GO 

Tau 

Tau-Aur , 
B22 

10.2 

5.9 

10 

13 

72 

DQ 

Tau 

Tau-Aur 

8.8 

6.4 

10 

13 

73 

Haro  6-37 

Tau-Aur 

8.8 

7.0 

10 

13 

74 

DR 

Tau 

Tau-Aur 

8.8 

5.7 

10 

13 

75 

DS 

Tau 

Tau-Aur 

8.8 

5.1 

10 

13 

76 

UY 

Aur 

Tau-Aur 

10.2 

2.7 

10 

13 

77 

GM 

Aur 

Tau-Aur 

8.8 

4.7 

10 

13 

78 

AB 

Aur 

Tau-Aur 

10.  2 

-2.  5 

10 

13 

79 

SU 

Aur 

Tau-Aur 

5.9 

3.4 

10 

208 

KV 

Mon 

NGC  2  2  64 

9.85 

6.1 

11 

210 

LL 

Mon 

NGC  2264 

9.85 

6.6 

11 

211 

PT 

Mon 

NGC  2264 

9.85 

7.6 

11 

212 

LM 

Mon 

NGC  2264 

9.85 

6.9 

11 

213 

NW 

Mon 

NGC  2264 

9.85 

5.7 

11 
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TABLE    III     (Continued) 


HRC 

STAR 

LOCATION 

(m-M) 

M 
PS 

REFERENCE    NOTE 

214 

LP  Mon 

NGC 

2264 

9.85 

6.7 

11 

215 

Lib  21 

NGC 

2264 

9.85 

4.7 

11 

216 

NX  Mon 

NGC 

2264 

9.85 

6.5 

11 

217 

W  84 

NGC 

2264 

9.85 

2.7 

11 

218 

LQ  Mon 

NGC 

2264 

9.85 

7.4 

11 

219 

LHa  2  5 

NGC 

2  2  64 

9.85 

3.3 

11 

220 

LR  Mon 

NGC 

2264 

9.85 

5.2 

11 

221 

LT  Mon 

NGC 

2264 

9.85 

5.9 

11 

222 

W  108 

NGC 

2264 

9.85 

2.5 

11 

223 

V419  Mon 

NGC 

2264 

9.85 

5.4 

11 

224 

V34  7  Mon 

NGC 

2264 

9.8  5 

7.6 

11 

225 

LU  Mon 

NGC 

2264 

9.85 

6.2 

11 

226 

10  Mon 

NGC 

2264 

9.85 

4.9 

11 

227 

IP  Mon 

NGC 

2264 

9.85 

4.6 

11 

228 

LY  Mon 

NGC 

2264 

9.85 

6.5 

11 

229 

LX  Mon 

NGC 

2264 

9.85 

5.7 

11 

230 

SS  Mon 

NGC 

2264 

9.85 

4.7 

11 

231 

V360  Mon 

NGC 

2264 

9.85 

4.1 

11 

232 

V4  2  6  Mon 

NGC 

2264 

9.85 

4.8 

11 

233 

V3  63  Mon 

NGC 

2264 

9.85 

6.2 

11 

234 

MM  Mon 

NGC 

2264 

9.85 

5.  2 

11 

235 

LHa  61 

NGC 

2  2  64 

9.85 

4.6 

11 

236 

V36  5  Mon 

NGC 

2264 

9.85 

5.  7 

11 

237 

V43  2  Mon 

NGC 

2264 

9.85 

6.1 

11 

238 

MO  Mon 

NGC 

2264 

9.85 

4.8 

11 

62 


TABLE  III  (Continued) 

HRC     STAR       LOCATION     (m-M)      M     REFERENCE    NOTE 
Pg     Pg 

23  9  OW  Mon 

24  0  OY  Mon 
241  MQ  Mon 
24  2  PY  Mon 

97  San  1 

98  HS  Ori 

99  HT  Ori 
100  V466  Ori 

102  P  1207 

103  VX  Ori 

104  VY  Ori 

105  San  2 

106  VZ  Ori 

107  V386  Ori 

108  P  1267 

109  P  1270 

110  Su  Ori 

111  WX  Ori 

113  P  1404 

114  EZ  Ori 

115  SW  Ori 

116  IU  Ori 

117  XX  Ori 

118  IX  Ori 

119  YY  Ori 
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NGC  2264 

9.85 

6.0 

11 

NGC  2264 

9.85 

5.0 

11 

NGC  2264 

9.85 

5.9 

11 

NGC  2264 

9.85 

4.7 

11 

Orion 

neb 

9.1 

3.9 

12 

Orion 

neb 

9.1 

6.6 

12 

Orion 

neb 

9.1 

5.3 

12 

Orion 

neb 

9.1 

4.6 

12 

Orion 

neb 

9.1 

5.6 

12 

Orion 

neb 

9.1 

6.7 

12 

Orion 

neb 

9.1 

6.8 

12 

Orion 

neb 

9.1 

4.9 

12 

Orion 

neb 

9.1 

5.5 

12 

Orion 

neb 

9.1 

7.5 

12 

Orion 

neb 

9.1 

5.4 

12 

Orion 

neb 

9.1 

3.8 

12 

Orion 

neb 

9.1 

6.7 

12 

Orion 

neb 

9.1 

6.1 

12 

Orion 

neb 

9.1 

3.1 

12 

Orion 

neb 

9.1 

3.3 

12 

Orion 

neb 

9.1 

5.6 

12 

Orion 

neb 

9.1 

1.1 

12 

Orion 

neb 

9.1 

6.7 

12 

Orion 

neb 

9.1 

5.5 

12 

Orion 

neb 

9.1 

5.2 

12 

TABLE  III  (Continued) 
HRC     STAR 


LOCATION 

(m-M) 

Pg 

M 
PS 

REFERENCE 

NOTE 

Orion 

neb 

9.1 

6.4 

12 

Orion 

neb 

9.1 

2.1 

12 

Orion 

neb 

9.1 

3.9 

12 

Orion 

neb 

9.1 

5.3 

12 

Orion 

neb 

9.1 

6.4 

12 

Orion 

neb 

9.1 

5.4 

12 

Orion 

neb 

9.1 

3.  5 

12 

Orion 

neb 

9.1 

6.1 

12 

Orion 

neb 

9.1 

7.1 

12 

Orion 

neb 

9.1 

6.2 

12 

Orion 

neb 

9.1 

5.7 

12 

Orion 

neb 

9.1 

5.8 

12 

Orion 

neb 

9.1 

2.8 

12 

Orion 

neb 

9.1 

4.0 

12 

Orion 

neb 

9.1 

4.9 

12 

Orion 

neb 

9.1 

6.0 

12 

Orion 

neb 

9.1 

5.7 

12 

Orion 

neb 

9.1 

4.8 

12 

Orion 

neb 

9.1 

5.6 

12 

Orion 

neb 

9.1 

6.3 

12 

Orion 

neb 

9.1 

5.6 

12 

Orion 

neb 

9.1 

5.8 

12 

Orion 

neb 

9.1 

7.3 

12 

Orion 

neb 

9.1 

7.0 

12 

Orion 

neb 

9.1 

4.6 

12 

120  YZ  Ori 

121  P  1649 

122  KM  Ori 

123  KN  Ori 

124  KP  Ori 

125  KR  Ori 

126  LL  Ori 

127  SY  Ori 

128  LN  Ori 

129  V356  Ori 

130  AA  Ori 

131  V486  Ori 

132  P  1817 

133  LX  Ori 

134  V488  Ori 

135  AB  Ori 

136  TT  Ori 

137  AD  Ori 

138  P  1931 

139  P  1946 

140  AG  Ori 

141  AI  Ori 

142  NS  Ori 

143  AL  Ori 

144  V360  Ori 
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TABLE  III  (Continued) 
HRC     STAR 


LOCATION 

(m-M) 

M 
Pg 

REFERENCE 

NOTE 

Orion 

neb 

9.1 

6.6 

12 

Orion 

neb 

9.1 

5.7 

12 

Orion 

neb 

9.1 

S.8 

12 

Orion 

neb 

9.1 

6.9 

12 

Orion 

neb 

9.1 

6.2 

12 

Orion 

neb 

9.1 

3.6 

12 

Orion 

neb 

9.1 

6.8 

12 

Orion 

neb 

9.1 

7.2 

12 

Orion 

neb 

9.1 

6.1 

12 

Orion 

neb 

9.1 

1.2 

12 

Orion 

neb 

9.1 

5.9 

12 

Orion 

neb 

9.1 

7.3 

12 

Orion 

neb 

9.1 

6.5 

12 

Orion 

neb 

9.1 

6.4 

12 

Orion 

neb 

9.1 

6.4 

12 

Orion 

neb 

9.1 

4.0 

12 

Orion 

neb 

9.1 

4.4 

12 

Orion 

neb 

9.1 

5.9 

12 

Orion 

neb 

9.1 

6.9 

12 

Orion 

neb 

9.1 

5.4 

12 

Orion 

neb 

9.1 

6.5 

12 

Orion 

neb 

9.1 

2.2 

12 

Orion 

neb 

9.1 

6.1 

12 

Orion 

neb 

9.1 

1.8 

12 

Orion 

neb 

9.1 

5.6 

12 

14  5  TW  Ori 

14  6  AM  Ori 

147  TV  Ori 

148  NY  Ori 

149  BO  Ori 

150  AN  Ori 

151  V573  Ori 

152  CE  Ori 

153  OT  Ori 

154  T  Ori 

155  AR  Ori 

156  V390  Ori 

157  AU  Ori 

158  V577  Ori 

159  AV  Ori 

160  PQ  Ori 

161  PU  Ori 

162  AZ  Ori 

163  Haro  4-125 

165  BD  Ori 

166  BC  Ori 

167  P  2441 

168  BE  Ori 

169  BF  Ori 
172  TX  Ori 
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TABLE  III  (Continued) 

HRC     STAR       LOCATION     (m-M)      M     REFERENCE    NOTE 
Pfi     Pfi 


173 

TY  Ori 

Orion 

neb 

9.1 

6. 

2 

12 

174 

Haro  7-1 

Orion 

neb 

9.1 

4. 

9 

12 

175 

Haro  7-5 

Orion 

neb 

9.1 

5. 

9 

12 

176 

Haro  4-2  55 

Orion 

neb 

9.1 

5. 

9 

12 

177 

V510  Ori 

Orion 

neb 

9.1 

5. 

6 

12 

178 

BH  Ori 

Orion 

neb 

9.1 

5. 

4 

12 

179 

Haro  7-4 

Orion 

neb 

9.1 

5. 

9 

12 

181 

DL  Ori 

Orion 

neb 

9.1 

4. 

4 

12 

182 

Haro  7-2 

Orion 

neb 

9.1 

3. 

9 

12 

188 

Lk  Ha  334 

Orion 

neb 

9.1 

5. 

,9 

12 

189 

Lk  Ha  33  5 

Orion 

neb 

9.1 

6. 

9 

12 

190 

Lk  Ha  336 

Orion 

neb 

9.1 

6. 

,9 

12 

191 

Lk  Ha  337 

Orion 

neb 

9.1 

4. 

.9 

12 

References : 

1.  Johnson,  1968 

2.  Johnson,  et  al.,  1961 

3.  Herbig,  1960 

4.  Buscombe,  1963 

5.  Herbig,  1971 

6.  Kholopov,  1959 

7.  Imhoff  and  iMendoza,  1974 

8.  Allen,  1973 

9.  Kuhi,  1964 

10.  Rydgren,  Strom,  and  Strom,  1975 

11.  Walker,  1956 

12.  Penston  and  Hunter,  1975 
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TABLE  III  (Continued) 

Notes : 

1.  ADiN  for  anonymous  dark  nebula. 

2.  Adopted  from  the  distance  modulus  for  FU  Ori  (HRC  186). 

3.  Assumed  AR  =  0.1. 

4.  Actually  4.0  =  My   not  M  . 

5.  Distance  modulus  adopted  from  the  average  modulus  calculated 
from  p  Oph  stars  discussed  in  Rydgren,  Strom,  and  Strom, 
1975. 

6.  Adopted  the  distance  modulus  for  V1057  Cyg  (HRC  300). 

7.  Adopted  the  distance  modulus  for  Cep  IV  as  that  of  Cep  III. 

8.  Adopted  from  average  distance  modulus  calculated  for  B209. 

9.  Adopted  from  average  distance  modulus  calculated  for  B7. 

10.  Adopted  from  average  distance  modulus  calculated  for  B19. 

11.  By  proximity  to  DI  Tau  (HRC  39),  adopted  the  same  distance 
modulus . 

12.  Adopted  from  average  distance  modulus  calculated  for  B18. 

13.  After  examination  of  Palomar  Sky  Survey  Prints  and  consider- 
ing the  absorptions  of  the  stars  in  Rydgren,  Strom,  and 
Strom,  1975,  the  star  was  assigned  to  one  of  three  groups 

on  the  basis  of  an  estimation  of  the  visible  absorption  in 
the  immediate  area.   An  average  distance  modulus  for  either 
light,  moderate  or  heavy  absorption  was  adopted. 

14.  Adopted  the  distance  modulus  based  on  the  average  of  the 
absorptions  in  Penston  and  Hunter,  1975. 
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9    10 
Absolute  photographic  magnitude  (M  ) 

Figure  27.   The  Luminosity  Functions,  <j)(MpP).   The  luminosity 

function  for  the  solar  neighborhood  is  adapted  from  Allen  [1973]. 

The  functions  for  the  other  four  plots  above  have  been  normalized 

such  that  the  number  of  stars  in  the  total  Orion  Population  is 

equal  to  the  number  found  from  integrating  (£(M  )  in  the  solar 

neighborhood  from  -2  <  M   <  10. 
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Figure  28.   H-R  diagram  showing  the  Hayashi 
track  for  a  one  solar  mass  protostar.   The 
bar  at  the  tracks  end  marks  the  stars  position 
on  the  main  sequence.   The  figure  is  adapted 
from  Clayton  [1968] . 
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Figure  29.   H-R  diagram  showing  the  Larson  track  for 
a  one  solar  mass  protostar.   The  solid  line  is  the 
zero- age  main  sequence'.   The  figure  is  adapted  from 
Strom,  Strom,  and  Grasdalen  [1975J . 
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